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ABSTRACT 


The PAMELA detector was launched on board of the Russian Resurs-DKl 
satellite on June 15, 2006. Data collected during the hrst four years have been 
used to search for large-scale anisotropies in the arrival directions of cosmic-ray 
positrons. The PAMELA experiment allows for a full sky investigation, with 
sensitivity to global anisotropies in any angular window of the celestial sphere. 
Data samples of positrons in the rigidity range 10 GV < R < 200 GV were 
analyzed. This article discusses the method and the results of the search for 
possible local sources through analysis of anisotropy in positron data compared to 
the proton background. The resulting distributions of arrival directions are found 
to be isotropic. Starting from the angular power spectrum, a dipole anisotropy 
upper limit S = 0.166 at 95% G.L. is determined. Additional search is carried 
out around the Sun. No evidence of an excess correlated with that direction was 
found. 


Subject headings: PAMELA, cosmic ray, anisotropy, positron 
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1. Introduction 

Measurements of cosmic-ray positrons address a number of questions in contemporary 
astrophysics, such as the nature and distribution of particle sources in our Galaxy, 
and the subsequent propagation of cosmic-rays through the Galaxy and the solar 
magnetosphere. Positrons are a natural component of the cosmic radiation, produced in the 
interaction between cosmic rays and the interstellar matter. Observations by the PAMELA 


experiment (Adrian! et al. 2009, 2010, 2013), conhrmed by other recent measurements 


(Ackermann et aL]|2012; Accardo et al.||2014), revealed a positron excess over the predictions 


of commonly used propagation models, e.g. (Moskalenko & Strong 2009 Delahaye et al 


2009). Either the standard paradigm of cosmic-ray propagation in the galaxy should be 


revised or additional sources of cosmic-ray should be introduced to explain the positron 
anomaly at high energy. There are two primary hypotheses for these new sources: 1) a 


source linked to particle physics, e.g. a dark matter decay or annihilation, e.g. (Girelli et al. 


2008). 2) a nearby astrophysical source, e.g. a pulsar or supernova remnant, e.g. (Hooper 


eFlI]|200^ 


In the hrst case, among the various dark matter models proposed. Bated et al. (2010), 


Schuster et al. (2010) suggested that a signihcant fraction of positrons detected at Earth 


could be produced by dark matter particles annihilating in the neighborhood of the Sun. In 
this case, the electron/positron pair emission would be strongly correlated with the center 
of the Sun and would produce visible effect in anisotropy study. 

In the second case, the localization of the astrophysical sources could results in 
anisotropy in the arrival direction of the cosmic-ray positrons. Several authors, e.g. 


(Ptuskin et al. 1995; Kobayashi et al. 2004 Bushing et al. 2008), estimated that the 


expected anisotropy in cosmic-ray positrons from supernova remnants and/or from pulsars 
should be of the order of percent or less. 
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Directional and timing data from PAMELA are exploited to search for positron 
anisotropies that can provide unique information on new sources of cosmic rays. Previous 
searches have been already carried out using the combined flux of electrons and positrons 


(Ackermann et ah 2012 Accardo et ah 2014 Campana et ah 2013) and results were 
compatible with an isotropic distribution of arrival directions of detected particles. 


In this paper the analysis is focused on positrons since they represent a cleaner sample 
of cosmic rays produced by possible sources respect to electrons. In most theoretical 
models, electrons and positrons from new sources are produced in pairs. However the e^ 
flux consists of only two components: e"*" from secondary production and from unknown 
sources, while the e“ flux contains, in addition and in a larger amount, also primary e“. 
Since in this energy range the cosmic-ray primary electrons represent the 90% of the total 
e“ flux, the search for anisotropies in the pure electron sample will be dominated by an 
isotropic background. 

In the following, an analysis method is presented to study the anisotropy of positrons 
using the back-tracing of particle trajectories, to identify arrival directions far from the 
Earth. An evaluation of systematic uncertainties is also provided. 


2. The PAMELA detector 


The PAMELA experiment has been operating on board of the satellite Resurs-DKl 
since June 2006. The design is optimized for precision studies of light particles and 
antiparticles in primary cosmic rays between a few tens of MeV and a few hundred of GeV. 


Since launch, PAMELA has collected ~ 8 ■ 10® events, of which ~ 5 ■ 10® electrons and 
positrons, whose analysis is described in several publications (Adriani et ah 2009, 2010| 


2011 ). 
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The PAMELA detector ( Picozza et al.||2007 ) is equipped with a magnetic spectrometer, 
comprising a permanent magnet hosting a tracking system. The tracking system consists 
of six double-sided microstrip silicon planes, which allow the determination of the particle 
charge and rigidity (momentum divided by charge) with high precision. An imaging 
electromagnetic calorimeter, made of 44 silicon planes interleaved with 22 plates of tungsten 
absorber, is mounted below the spectrometer to accurately perform particle identihcation 
and to measure the energy of electrons and protons. A Time of Flight (ToF) system made 
of three planes of double layers of plastic scintillator permits measurements of particle 
velocity and energy loss. It also provides the main trigger for the experiment. Another 
layer of scintillator and a neutron detector, placed below the calorimeter, give additional 
information about the shower extent and improves lepton/hadron discrimination. Particles 
in the PAMELA acceptance due to scattering or interactions are rejected during off-line 
analysis of the anti-coincidence system signals. 


This apparatus permits electrons and positrons to be separated from the proton 
background over the rigidity range 10 GV < R < 200 GV and to measure their incoming 
direction with an accuracy of about two degrees. For each particle crossing the detector the 
arrival direction from space is reconstructed using the trajectory inside the instrument and 
the satellite position. The satellite orbit, 70° inclination and 350-610 km altitude, allows 
PAMELA to perform a very detailed measurement of the cosmic radiation in different 
regions of Earth’s magnetosphere. Back tracing of detected particles in the geomagnetic 
held is performed in order to obtain their initial spatial distribution, outside of the Earth 
magnetosphere. 
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3. Analysis Method 


This data analysis aims to identify the presence of any large-scale pattern in the 
distribntion of arrival directions of cosmic-ray positrons detected by PAMELA. The entire 
data-set of PAMELA np to Jannary 2010 was analyzed, constrncting sky maps for positrons 
in Galactic coordinates. In order to perform a search for anisotropies in these data, the 
experimental particle flnx, measnred by the instrnment for different directions, is compared 
with backgronnd (or coverage) maps derived from an isotropic particle flnx. 


Coverage maps can be determined throngh an accnrate simnlation of particles coming 
from any direction of the sky, bnt this wonld reqnire a very precise knowledge of the 
instrnment exposnre, defined as the projection of the acceptance in each sky direction 
integrated over the fnll live time of the detector. Since the resnlt depends on the position 
and orientation of the spacecraft as a function of time, small inaccnracies can lead to the 
creation of fake anisotropy signals. To overcome this problem, the backgronnd map for 
isotropic expectation was constrncted nsing PAMELA data themselves. The flnx of protons 


is highly isotropic, as previonsly shown with PAMELA data (Giaccari et ah 2013). A set of 
proton events was selected as described in Sec. |^for the same period of time as nsed for the 
positron data. The instrnment exposnre, dead times and other detector effects are therefore 
inclnded when calcnlating the ratio of positrons over protons. The comparison between 
signal and backgronnd maps was performed with two independent and complementary 
techniqnes: a statistical significance test introdnced by Li and Ma ( |Li et ah 1983) and a 
spherical harmonic analysis. 










3.1. Back tracing 


In order to identify particle arrival directions and investigate possible anisotropies, 
it is necessary to account for the effect of the geomagnetic held on particle propagation. 
Consequently, trajectories of all selected events were reconstructed in the Earth’s 


magnetosphere by means of a tracing program (Bruno et ah 2014) based on numerical 


integration methods (Smart et ah 2000). The IGRF-11 (Finlay et ah 2010) model and the 


TS05 (Tsyganenko et al. 2005) model were used for the description of the geomagnetic 
held: the former employs a global spherical harmonic implementation of the main magnetic 
held; the latter provides a dynamical (hve-minute resolution), semi-empirical best-ht 
representation of the external geomagnetic held sources in the inner magnetosphere, based 
on recent satellite measurements. Using satellite ephemeris (position, orientation, time 
information) and the particle rigidity and direction measured by the tracking system, 
trajectories were back-propagated to the model magnetosphere boundaries, allowing 
corresponding asymptotic arrival directions (i.e. the directions of approach before they 
encountered the geomagnetic held) to be determined. 


3.2. Sky Pixelation 


To determine the arrival direction distribution over the whole sky, the Healpix 
framework was used. Healpix provides a visualization of the sky map in a 2D-sphere 


(Gorski et al. 2005). This tool pixelises (i.e. a subdivides) the 2D-sphere, where each pixel 


covers the same surface area as every other pixel and is regularly distributed on lines of 
constant latitudes. In the Healpix scheme the pixel size cannot be chosen arbitrarily, the 
total number of pixels is given by 12 x where Ugide dehnes the map resolution and 
can take only a value that is power of two. A grid with 3072 pixels with an angular pixel 
extension of about ~ 7° was adopted. All maps are given in Galactic coordinates. 
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3.3. Significance Maps 


The statistical significance of the signal over an isotropic backgronnd in each direction 


of the sky was derived following the techniqne described in the work of Li and Ma (Li et ah 


1983). The signihcance formula is dehned as 


5 = ±y2<^ Ar„„ln 


iVoffln 


1 T n 


a 




(1 + a) 


Non + Nofl 

Non 


+ 


Non + iV. 


off 


1/2 


( 1 ) 


where Non and iVofr are respectively the observed and the expected number of events in a 
certain angular window of the sky, after the proton and positron maps are normalized to 
the same number of events. 


In Eq.([^ a is the ratio of the exposure in the on-source region divided by the exposure 
in the off-source region; in this case the on-source and off-source regions coincide so a = 1 . 
The signihcance is dehned as positive if there is an excess otherwise as negative in the 
presence of a dehcit of events. If no signal is measured, the distribution of the signihcance S 
is a normal function. Results are reported bin per bin, as a number of standard deviations. 

This technique, developed in 7 -astronomy to search for point sources, can be easily 
arranged to search large anisotropies integrating the signal and the background maps over 
various angular scales in order to increase the sensitivity up to a large-scale pattern. In this 
survey diherent angular radii for the integration, from 10° to 90°, were used. The center 
of each bin of the sky pixelization is considered and the content of the bin is given by the 
cumulative number of events falling within a given angular distance from the center of that 
bin. Both signal and background maps are smoothed following this process. In this way the 
obtained integrated maps have correlated neighboring pixels and none of the information at 
a given angular scale is lost. 
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The statistical Li and Ma significance test is subsequently applied to study possible 
deviations from isotropy. In the case of no strong anisotropies, each significance histogram 
approximates a Gaussian distribution. The significance distributions becomes narrower 
with increasing integration radius - this is expected as in the integrated maps the bins are 
strongly correlated and identical events are included in the calculation. 


3.4. Spherical Harmonic Analysis 


The orbital parameters of the Resurs-DKl satellite provide a relatively uniform 
exposure over the full celestial sphere. This permits the study of the angular power 
spectrum of arrival directions of the events, i.e. the cosmic-ray intensity over the celestial 
sphere can be expanded in spherical harmonics. The power spectrum analysis is a powerful 
tool to study the anisotropy patterns and provides information on the angular scale of the 
anisotropy in the map. 


The positron (or electron) relative intensity map I, defined as the relative deviation 
of the observed number of events N from the expected number of events < iV > in each 
angular bin of the sky defined by the Galactic coordinates {glon, glat) is considered. 


I{glon, glat) 


N{glon, glat)— < N{glon, glat) > 
< N{glon, glat) > 


( 2 ) 


The relative map I is expanded in the basis of spherical harmonics functions Yim 

oo m=l 

I {glon, glat) = EE aimYim{glon,glat). (3) 

/=0 m=—l 

The full anisotropy information is encoded into the set of spherical harmonic coefficients, 
aim, which can be used to calculate the angular power spectrum 

m=l 

Ci') = ^ E “L. 

m=—l 


( 4 ) 
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The amplitudes of the power spectrum at some multipole order are associated to the 
presence of structures in the distribution of the arrival directions in the sky. Non-zero 
amplitudes in the multipole moments, Ci, arise from fluctuations in the particle flux on an 
angular scale near 180// degrees. 


3.4-1- Dipole upper limit 


Although all the spherical harmonic moments are important to characterize the 
anisotropy patterns at any angular scale, the dipole moment is of particular interest 


(Bushing et ah 2008). In fact, if there is a marked directionality, dipole anisotropy is 


expected as could be the case for a single source dominating the positron flux. In this case 
the overall intensity at an angular distance 9 from the maximum of the dipole anisotropy 
can be written as I(6')=lo+Iicos(6'), where Iq represents the isotropic signal and Ii represents 
the maximum anisotropy due to the dipole. The degree of anisotropy can be expressed as 
the fraction 5 = Ii/Iq. 


By applying Eq. the intensity of the dipole anisotropy becomes 

I{e)- < 1(9) > _ I{9) -lo _ hcos{9) 


m = 


< m > 


(5) 


Since the spherical harmonic Yi{9,(j)) depends from 9, it is possible to derive that 


a? = ^ and, consequently. 


_ //i\247r 


( 6 ) 


Starting from Eq. the correlation between the dipole amplitude and the power spectrum 
coefficient with / = 1 is 

frD 

(7) 


V dvr 


Considering Ci as the outcome of the calculated angular power spectrum over the signal 
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map for /=1, for the observed spectrum 

V dvr 

Since Ci follows a distribution with three degrees of freedom centered on Ci, the 
probability to measure Ci given the true Ci is 


( 8 ) 




3^3 Cl 


exp 


3Ci 


V^Ci V C'l " V 2Ci 

From Eq. and it is possible to derive the probability for 5 


( 9 ) 


as 


( 10 ) 


The upper limit on 6 can be calculated using a frequentist approach where the limit 
corresponds to the solution of the integral 

rSm 

/ P{6\6) d5 = 1-C.L. 


( 11 ) 


for a given conhdence level (C.L.) (Neyman 


1937 


Garthwaite et ah 


1995), where (5n 


represents the dipole anisotropy corresponding to the measurement of the dipole power 
spectrum Ci^rn- 


4. Data Selection 

The sample of positrons selected for this analysis concerns to the period from the 
satellite launch to the end of the solar minimum phase (June 2006 - January 2010). This 


is the same data set as used for previous publications (Adriani et al. 2010). The selection 


criteria are, however, different since a statistical subtraction of proton contamination in 
the positron sample does not allow event-by-event positron identihcation as required to 
determine the incoming direction and rigidity of each event. The positron sample was 
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instead selected using a cut based analysis similar to the one developed for the electron 


flux determination (Adriani et ah 2011). The same track and event quality selection 
were used; the hadron contamination was reduced to a negligible amount by using a 
stronger calorimeter selection (with about 80% selection efficiency) combined with activity 
requirements on the neutron counter. A very clean sample of positrons was obtained from 
flight data in the rigidity range from 10 to 200 GV. A total number of 1489 positrons with 
rigidity R > 10 GV was selected for this work. With the same selection a sample of 20673 
electrons and positrons was selected for the anisotropy search in the direction of the Sun. 
The quality of the sample was verified by comparing the resulting positron fraction to the 
published one. In the same period, June 2006 to January 2010, the calorimeter was used 
to select hadronic showers initiated by positively charged particles in the same rigidity 
range and with the same requirements on the event and track quality. In this way a proton 
sample of about 4.5 ■ 10^ events, which preserved the angular distribution of particles along 
PAMELA orbit, was obtained. Finally, a refined analysis on the detector orientation was 
carried out by comparing the pointing of the instrument determined by the satellite orbital 
information with the pointing obtained using data from the Earth remote sensing camera 
that is installed on board the Resurs-DKl satellite. This refined analysis minimized the 
uncertainty on the measured absolute incoming trajectory for all the selected particles 
which is used in the back-tracing and in the anisotropy maps. 


5. Results 

Fig. HI shows the maps obtained for events with rigidity greater than 10 GV, for 
selected positrons and protons, taking into account the geomagnetic field. As the PAMELA 
experiment has a shorter exposure to the poles, the sky is not observed uniformly, leading 





14 



(a) Positrons (b) Protons 

Fig. 1.— Event maps for positrons (a) and protons (b) for 10 GV < R < 200 GV, taking 
into account the geomagnetic field effect. 


to a particular shape in the observed maps. The same overall exposure is achieved for 
both types of particle. Both the signal and the background maps were binned taking into 
account the detector angular resolution. The color scale represents the number of detected 
events for each bin. If a signal can be found in the data, it could spread over multiple 
adjacent bins, making the anisotropy measurement more difficult. Hence, to highlight a 
possible anisotropy, the bin size must be similar to the angular scale of the anisotropy itself. 
To increase sensitivity it is possible to merge the content of neighbouring bins. The content 
of each integrated bin is equal to the total number of events belonging to the bins covered 
by a circular region of a given radius. The radius chosen for the integration represents the 
angular scale at which the anisotropy is expected. If the integration radius is too small the 
signal is distributed among adjacent bins, while if it is too large too much background will 
be integrated. The events and the coverage maps were integrated over the following angular 
scales: 10°, 30°, 60°, 90°. The corresponding results are reported in Fig. [^for positrons 
and in Fig. [^for protons. The color scale now represents the integrated number of events. 
The exposure distribution varies depending on the integration radius as expected. To 


estimate the statistical signihcance of any excess or dehcit the method detailed in Sec. 3.3 


was used. In this method, likelihood functions are applied both to the null hypothesis and 
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Fig. 2.— Positron maps for 10 GV < R < 200 GV, over the following angular scales: 10°, 
30°, 60°, 90°. 

to the signal one; the ratio between these two functions represents the signihcance of the 
excess. Signihcance maps, constructed by comparing signal and corresponding background 
maps for the different integration radii, are shown in Fig. The signihcance signals are 
always less than four sigma, indicating that a point-like source cannot be distinguished from 
background. In the case of statistical huctuations for an isotropic sky the distribution of the 
signihcance will exhibit a Gaussian distribution with zero mean. In Fig. the obtained 
signihcance distributions are shown together with a htted Gaussian function; no signihcant 
deviation from isotropy is observed. It can be seen that the signihcance distribution 
becomes narrower when increasing the integration radius since with larger integration 
radius, the bins are strongly correlated, reducing the variance of the distribution. For the 
study of the angular power spectrum of the arrival direction distribution, the cosmic-ray 
intensity was expanded in spherical harmonics, using the anafast code provided by Healpix 























































16 






Fig. 3.— Proton maps for 10 GV < R < 200 GV, over the following angular scales: 10°, 
30°, 60°, 90°. 


software (Gorski et al. 2005). The power spectrum in modes from I = 1 (dipole) up to 


I = 20 was studied. Due to statistical and systematic uncertainties higher order modes are 
irrelevant for this search. In Fig. |^the angular power spectra G(l) as a function of the 
multi-pole 1 are shown. The band dehned in the hgures by the two dotted lines shows the 
3 ct region of the expected power spectrum from an isotropic sky. The power spectra are 
compatible with the prediction of the isotropic sky and no signihcant deviation from the 
isotropy was observed. The value for the measured dipole amplitude (5i, as expressed in Sec. 
3.4.11, is calculated with the value of coefficient Ci reported in Fig. Solving the integral 


in Eq. ([IT|) for a conhdence level at 95%, the limit on the dipole anisotropy parameter is 5 


0.166. 
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Fig. 4.— Significance maps for 10 GV < R < 200 GV, over the following angnlar scales: 
10°, 30°, 60°, 90°. 


6. Systematic uncertainties estimation 

Systematic uncertainties were estimated with the following methodology. Variables 
that have an effect on the anisotropy measurement were independently and artihcially 
varied in the selected data samples to understand the effects of the introduced biases 
on the anisotropy measurement. For each of the new biased samples of positrons and 
protons the analysis was performed as explained in the previous sections; the resulting 
Li and Ma distribution and angular power spectrum were compared to the un-biased one 
and the difference was taken as an estimation of systematic uncertainties for the given 
bias. Assuming no correlation between the systematic effects (biases), the total resulting 
systematic uncertainty was obtained by quadratically summing the single components. 


Three main sources of systematic errors were identihed: 











































Number of pixels Number of pixels 


~ 18 - 




(a) (b) 




(c) (d) 

Fig. 5.— Histograms for the significance maps for 10 GV < R < 200 GV, over the following 
angular scales: 10° (a), 30° (b), 60° (c), 90° (d). The grey area represents the systematics. 


1. particle direction reconstruction. Experimentally, the particle direction reconstruction 
in celestial coordinates results from the combination of the particle direction in the 
PAMELA reference frame measured by the tracking system and the pointing direction 
of the satellite (hence of PAMELA which is solidly attached to the Resurs-DKl). The 
inclination information of the satellite is provided to PAMELA by the satellite every 
1.5-180 s. The higher rate is used during the satellite movements that are performed 
to observe the Earth with the optical camera mounted on the satellite. Hence, satellite 
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Fig. 6.— Angular power spectra C(l) as a function of the multi-pole 1 for the positron signal 
over the proton background. The grey area represents the systematic effects calculated as 
described in Sec. 1 

inclination information derives from an interpolation of points as function of time 
and it has a resolution usually better than about two degrees. To estimate the effects 
of a possible systematic in the particle direction reconstruction on the anisotropy 
measurements, the inclination measurement was changed for each event by choosing 
random values according to a Gaussian distribution, with a variance equal to the 
experimental angular resolution. 

2. Particle rigidity determination. Particle rigidity is measured by the tracker, by 
determining the deflection of charged particles in the magnetic held. The maximum 
detectable rigidity (MDR) in PAMELA is of about 1.2 TV, resembling the accuracy 
in the cluster position determination of about 3 /xm in the bending view and the 
spectrometer magnetic held of about 0.4 Tesla. Residual coherent misalignment of the 
tracker silicon sensors after detector calibration could bring a systematic shift in the 
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rigidity measurement of about 15% at the MDR corresponding to a maximum shift in 
the deflection measurement of about 0.0001 GV“^. A systematic error in the rigidity 
measurement could have an impact in the anisotropy measurement, since the applied 
back tracing procedure relies on the particle rigidity to perform the integration of 
motion in the Earth’s magnetic held. To estimate the effects of these systematic 
uncertainties, a dehection uncertainty ±10“"^ GV“^ was added and subtracted to the 
dehection measurement of selected particles. 

3. Event time determination. Absolute UTG time is saved by the PAMELA GPU with a 
precision of one millisecond. The absolute time of the PAMELA GPU is synchronized 
with the satellite with a precision of one second at least once per orbit. Different clock 
frequencies between the two clocks and the absolute time at ground are negligible. 
Time determination is used to calculate the satellite position in latitude, longitude 
and altitude given the orbital elements. Hence, the uncertainty in time corresponds 
to an uncertainty in the orbital positioning and inhuences the back tracing procedure 
and the anisotropy estimation. Systematic effects were estimated by biasing the 
sample adding and subtracting one second to the measured time on an event-by-event 
basis. 

The resulting estimated systematic uncertainty on the anisotropy measurement is of 
the order of about 12% and shown as grey bands in hgures and 


7. Search in the direction of the Sun 

An excess of electrons and positrons was also searched for in the Sun direction. The 
reference frame used in this analysis is the ecliptic coordinate system centered on the Sun 


and referred to the J2000 epoch (Murray 


1989). 
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The number of electron and positron events within annuli centered on the Sun in the 
range 0°-90°, in 5° steps, is calculated. The results are reported in Fig. [^as a function of 


S 700- 

O I 

I 600 E- 

E : 

i 500^ 

400 ^ 

300 ^ 

200 ^ 

100 - ♦ 


• Observed 
I I Expected 


60 80 100 
Angular distance to Sun (°) 


Fig. 7.— Total number of electrons and positrons with 10 GV < R < 200 GV as a function of 
the angular distance from the Sun. The grey boxes correspond to the 3a fluctuation respect 
to an isotropic expectation. 


the angular distance from the Sun direction and are compared to the expectations from 
an isotropic sky. The grey boxes represent the 3cr bounds for an isotropic sky calculated 
with an elliptical reference frame on subsets of coverage map. Data are consistent with 
the isotropic expectation within a 3cr interval. Similar results were obtained by Fermi LAT 


collaboration (Ajello et al. 2011) without selecting electrons and positrons separately. 


8. Discussion 


The detection of an anisotropy signal might be used to distinguish between models that 


could explain the cosmic-ray positron excess, e.g. (Adriani et al. 2009), over the predictions 


of commonly used propagation models. A variety of models, invoking either dark matter or 


astrophysical objects as positron sources, e.g. see review papers (Boezio et al. 2009 Serpico 


2012), have been put forward to explain this excess. As a result of the different localization 
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of such sources an anisotropy in the arrival particle distribution cannot be excluded. 


Various authors, e.g. (Ptuskin et ah 1995 Ackermann et ah 2010; Di Bernardo et ah 2011 


Kobayashi et ah 2004| Bushing et al. 2008), studied the expected anisotropy in cosmic-ray 


positrons from supernova remnants and/or from pulsars. For example, Ackermann et al. 


(2010), Di Bernardo et al. (2011) estimated an anisotropy of the order of about 1% for 


100 GeV positrons accelerated by Vela supernova remnant, while Bushing et al. (2008) 
estimated an anisotropy of about 3% and about 0.3% for positrons with energies between 
10 and 100 GeV accelerated by Monogem and Geminga pulsars, respectively. In the case a 


dark matter origin, Ackermann et al. (2010) estimated an anisotropy of the order of per-mil. 


A blind search for positron anisotropies was performed on PAMELA data, taking into 
account the effects of the Earth’s geomagnetic held. The analysis of PAMELA data was 
carried out by using two techniques: the signihcance maps, which allowed the study of the 
directions of possible over-densities at any angular scale, not just the dipole contribution, 
and the angular power spectrum, which is a method capable of determining the anisotropy 
magnitude. Results were consistent with isotropy at all angular scales considered. A limit 
on the dipole anisotropy parameter, S = 0.166 at 95% conhdence level, was obtained. This 
limit is close but higher than the expected anisotropy for astrophysical sources and about an 
order of magnitude higher than in dark matter models that, therefore, cannot be excluded. 
Because of the limited statistics of the positron data set it was not possible to investigate 
the energy dependence of the dipole amplitude. 


A few models (Batell et al. 2010 Schuster et al. 2010) predict that large astrophysical 
bodies within the Solar System may signihcantly contribute to cosmic-ray positrons 
capturing dark matter particles that may annihilate hrst into metastable mediators. If the 
mean free path of these mediators is in excess of the solar (or planetary) radius, they can 
escape and annihilate again, producing gamma rays and charged particles such as positrons 

















































23 


arriving from directions correlated with the centers of the astrophysical source. The analysis 
of PAMELA all electron data in the direction of the Sun are consistent with the isotropic 
expectation at a 99.7% conhdence level. 


Results of the analysis reported in this paper are in agreement with those published 
by AMS and Fermi LAT. It should be noted, however, that a special care is needed when 


comparing the results of these three different experiments. AMS results (Accardo et ah 


2014) refer to positrons but reporting only the dipole upper limit, while in this paper a 


multipole analysis is presented. Fermi LAT results (Ackermann et ah 2012) refer to the 


sum of electrons and positrons. In the energy range between tens of GeV and hundreds of 
GeV, electrons from standard cosmic-rays sources are the dominant component. Hence, the 
study of the pure positron signal is more sensitive, in this energy range, to anisotropies. 
Finally, in this analysis the effects of the Earth’s magnetic held, which could affect a weak 
anisotropy, were taken into account, while both the AMS and Fermi LAT results refer to 
the anisotropy as measured in a low Earth orbit. 

It should be noted that the local (a few tens of parsec) turbulent realization of the 


interstellar magnetic held (Giacinti et ah 2012) as well as magnetic reconnection in the 


heliotail ( Lazarian et al.||2010 ) may have an ehect on the cosmic-ray arrival directions. Such 
processes have been proposed to explain the anisotopy observed mostly in the multi-TeV 


energy region (Nagashima et ah 1998; Amenomori et ah 2006 Guillian et ah 2007; Abdo et 


al. 2009 Aglietta et ah 2009 Abbasi et ah 2010 luppa et ah 2013). However, there are no 


quantitative studies on the ehects of these processes on the arrival directions of tens of GeV 
positrons measured at Earth. 
























































9. Conclusions 


The arrival direction of cosmic-ray positrons in the rigidity range from 10 GV to 
200 GV for about 10^ positrons collected in the hrst 4 years of operation of PAMELA 
was studied. The effects of the Earth’s geomagnetic held were accounted for. Results are 
consistent with isotropy at all angular scales considered. Moreover, no enhancement of the 
all electron hux was observed in the direction towards the Sun. 
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